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Nonperturbatively Improved Hadron Spectroscopy Near the Continuum 
Limit 
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We report the results of our quenched lattice simulations of the Wilson action with a nonperturbatively de- 
termined clover term at /3 = 6.2 and compare them with those of the standard Wilson action at the same /3 
value. 



Simulations in full QCD are extremely slow, 
and are still carried out far from the continuum 
limit. A possible solution for this problem is the 
use of improved actions, for which the simula- 
tion gives the continuum limit at higher values of 
the cutoff a, reducing computational costs. An 
important step in the improvement program was 
taken recently by the Alpha collaboration, with 
the nonperturbative determination of the coeffi- 
cient cc;]y in the clover action . 

We present here the results of our study of 
hadron spectroscopy using the clover-improved 
Wilson action (sec for notation) 
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The nonperturbative expression of csw is given in 
terms of the coupling /3 in |^ . We consider lattice 
volume 24P x 48 and coupling /3 — 6.2. We thus 
take Csw = 1-61375065. We also consider the 
unimproved case {csw = 0) for comparison. For 
the improved case we consider the following val- 
ues for the hopping parameter k: 0.1240, 0.1275, 
0.1310, 0.1340, 0.1345, 0.1350, 0.1352. For the 
csw = case we consider: 0.1350, 0.1400, 0.1450, 
0.1506, 0.1510, 0.1517, 0.1526. The simulations 
were carried out on the APE 100 computer at Tor 
Vergata. 

Our statistics come from 50 quenched gauge 
configurations, generated by a hybrid overrelax- 
ation algorithm with each update corresponding 
to a heat-bath sweep followed by 3 overrelaxation 



sweeps. The configurations are separated by 1000 
updates. 

The inversion of the fermion matrix is per- 
formed using the stabilized biconjugate gradient 
algorithm 1^. We restart the inversion from the 
current solution every 100 iterations, in order to 
reduce accumulation of roundoff errors Q. We 
employ point-like sources. 

We sum fermion propagators over the space 
directions x,y,z for sites within blocks of side 
3, and then store the result. We then form 
hadron correlations from these "packed" propaga- 
tors. This procedure becomes exact in the limit of 
an infinite number of configurations, since it cor- 
responds to including gauge-noninvariant terms 
in the computation of hadron correlations (in our 
case we have checked that the errors thus intro- 
duced are negligible). This corresponds to gain- 
ing a factor 3'^ in storage, and has enabled us to 
have all quark propagators stored simultaneously. 
In this way we can form hadron correlations from 
nondegenerate combinations of quark flavors, as 
done in Q. 

Within the improvement program one can de- 
termine the critical value Hc using the unrenor- 
malized current quark mass or "Ward identity" 
mass, defined through the improved PCAC rela- 
tion as |l| 
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Linear extrapolation to the limit of mwi — 
provides a much more stable fit for the dctermi- 
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nation of than the conventional fit of pseudo- 
scalar meson masses to the limit of zero pion 
mass. Our results for Kc are given below. 



improved case: 

from mwi 
from 
cgw = case 



Mis 



(from Ml 



: Kc = 0.135802(6) 
: Kc = 0.135861(19) 

0) = 0.153307(19) 



Once the value for has been determined, 
one can obtain chiral extrapolations from plots 
of hadron masses as linear functions of the bare 
quark mass, defined as mg^K) = (1/k — 1/kc)/2, 
or as linear functions of the measured squared 
pseudoscalar mass Mpg(K). In the improved case 
one can define the modified or "improved" bare 
quark mass ||l| by rriq (k) = mq{K) [1 + bm mq{K)] 
(Note that fhq is the renormalized mass with 
Zjri = 1 ) The improvement coefficient has 
been determined nonperturbatively 6m — 
—0.60 . Note that since we have nondegenerate 
flavor combinations we use averages of the masses 
defined above, e.g. for a meson corresponding to 
the flavors ki and K2 we define mq{Ki, K2) = 

[mq{Kl) + mg(K2)]/2. 
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Table 1 


Values of in MeV. 






Exp. 


IMPROVED 


csw = 


Jk' 


0.487 


0.40(14) 


0.40(9) 




0.557 


0.44(16) 


0.45(10) 



Table 2 : Values of the quantity J and comparison 
with experiment. 

In Table we give our values of the inverse lat- 
tice spacing coming from ratios of the vector 
meson mass My over its experimental value at 
different fiavor combinations. We determine the 
strange-quark mass in lattice units from two 
methods: a fit of {Mps/M^Y (method I) and 
a fit of [Mps /MvY (method II) as functions of 
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Figure 1. Octet mass in terms of (i)A/pg (ii)?^^. 

the quark mass, interpolating to the experimen- 
tal values of these ratios. (Here My denotes the 
chiral extrapolation of My.) For the improved 
case we use fhq as the bare quark mass, while in 
the unimproved case we use mq. In Table | are 
the values of the quantity J = My {dMy /dM^g) 
obtained from {Mr/Mk-Y vs. M^g. 

All our baryon mass values for the improved 
case come from linear fits of the masses as func- 
tions of the improved quark mass m^, except 
for the S — A mass splitting, where we used a 
quadratic fit. (Similarly for the csw = case, 
using as variable mq.) 

Note that we have defined the quantity Moct = 



3 



{Mm + M\)/2 (this combination was chosen so 
that the resulting mass is flavor-symmetric). For 
this quantity the advantage of using the improved 
quark mass is clearest, see Fig. |^. (A similar im- 
provement is observed also with respect to plots 
using the unimproved bare quark mass ruq.) 

We show an APE plot in Fig. || for the im- 
proved case. Note that we divide at each point 
by Mv interpolated to the strange-quark mass, 
corresponding to M$ (i.e. a constant value) . Ex- 
perimental points in these figures correspond to 
Mtv, My: a-nd M-^ and appropriate meson masses. 
In Table m we give our mass values in MeV. Table 
H contains data in lattice units for the improved 
case (diagonal flavor combinations only). 



Particle 


Exp. value cgvy = 


IMPROVED 


Mat 


939 


977(115) 


953(112) 


Ma 


1115.7 


1142(113) 


1127(115) 


Ma 


1232 


1382(172) 


1307(160) 


Moct 


1152.55 


1138(113) 


1141(113) 


Ma-n 


293 


381(90) 


366(86) 




232.45 


320(70) 


295(66) 




73.7 


50(39) 


59(25) 


Table 3 : Baryon masses in MeV 
with experiment. 


and comparison 


K 


Mps" 


My 


Mat 


0.1240 


1.15951 


1.1044(18) 


1.7046(54) 


0.1275 


0.73101 


0.8922(19) 


1.3821(50) 


0.1310 


0.367794 


0.6595(22) 


1.0220(55) 


0.1340 


0.117867 


0.4379(42) 


0.6588(91) 


0.1345 


0.0833452 


0.3995(59) 


0.5929(100) 


0.1350 


0.0517096 


0.3660(133) 


0.5169(141) 


0.1352 


0.0397285 


0.3543(225) 


0.4859(177) 



Table 4 : Masses in lattice units. 



With our present statistics the main effect of 
the nonperturbative improvement observed on 
hadron masses is a smaller spread in the value 
of the lattice spacing extracted from mesons with 
and without strange quarks. We have also shown 
how the use of the improved quark mass turns the 
rough behavior of the dependence of the octet 
baryon mass upon quark masses into a smooth 
one. 

We are currently increasing our statistics 0, 
and we plan to extend our simulations to the case 
with bermions ||], an approximate method that 
allows for an estimate of dynamic quark effects. 
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Figure 2. APE plot for the nucleon mass. 
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